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We investigate the presence of multipartite entanglement in macroscopic spin chains. We discuss 
the Heisenberg and the XY model and derive bounds on the internal energy for systems without 
multipartite entanglement. Based on this we show that in thermal equilibrium the above mentioned 
spin systems contain genuine multipartite entanglement, even at finite modest temperatures. 
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I. INTRODUCTION 

The study of entanglement in condensed matter sys- 
tems has become a major line of research in quantum 
information science. In these studies, different aspects 
of entanglement have been investigated. First, by inves- 
tigating the entanglement in the reduced state of two 
qubits in a spin chain, it has been shown that entangle- 
ment theory can help to understand the physical proper- 
ties of systems undergoing quantum phase transitions . 
A slightly different approach studies the entanglement 
between two qubits which can be generated by local op- 
erations on the remaining qubits, the so-called localizable 
entanglement 01 • A third field of study investigates the 
entanglement between a block and the remaining qubits 
in the ground state. The entanglement can be quantified 
by the entropy of the reduced state and then the question 
arises, whether the entanglement scales like the surface 
of the block 0. Furthermore, the insights from entan- 
glement theory made it possible to understand the limits 
of known simulation techniques (like the density matrix 
renormalization group) in condensed matter physics as 
well as the design of new techniques which are superior 
to the ones known before 4]. Finally, it has been shown 
that macroscopic properties of solids can be related to en- 
tanglement properties of microscopic degrees of freedom, 
allowing for a detection of entanglement in macroscopic 
objects by observing macroscopic observables |5|, |y,l3- 

In all of the research lines specified above, entangle- 
ment was studied as a bipartite phenomenon; the inves- 
tigation of multipartite entanglement in spin chains has 
only a limited literature 0, I n entanglement theory, 
however, multipartite entanglement has been intensively 
studied and it is known that quantum correlations in the 
multipartite setting have a much richer structure than in 
the bipartite setting jToL H3- The main results known 
are for three or four particles, results for large ensembles 
of particles are still rare. Nevertheless, studying multi- 
partite entanglement in large systems can already benefit 
from the results for small systems in two respects. First, 
on a more local scale, one can pick up sub-ensembles of 
three or four particles and investigate the multipartite 



entanglement of the corresponding reduced state. Sec- 
ond, on a more global scale, one can ask what types of 
multipartite entanglement are necessary to form a given 
state of the total system. 



In this paper we investigate multipartite entanglement 
in macroscopic spin systems. The tool we use is simi- 
lar to the one introduced in Ref. 0: We derive bounds 
for the internal energy U = (H) which have to hold 
for states without multipartite entanglement, violation 
of these bounds implies the presence of multipartite en- 
tanglement in the system. We investigate in detail two 
one-dimensional spin models: the Heisenberg model and 
the isotropic XY model. Our results lead to the insight 
that for these models even at modest temperatures the 
correlations cannot be explained without assuming the 
presence of genuine multipartite entanglement. However, 
we would like to stress that in our approach we do not 
assume that the spin system is in the ground state or in a 
thermal state, our theorems can also be applied to states 
out of equilibrium. 



Our paper is organized as follows: In Section II we 
introduce the basic definitions of multipartite entangle- 
ment we want to apply to spin systems. We first re- 
call the notion n-separability and genuine multipartite 
entanglement. Then, we also introduce the notion of 
/c-producibility, which is well suited for the investiga- 
tion of multipartite entanglement in macroscopic sys- 
tems. In Section III we apply these terms to the anti- 
ferromagnetic Heisenberg chain. We calculate energy 
bounds below which either reduced states are genuine 
multipartite entangled or the total state requires multi- 
partite entanglement for its creation. This shows that 
already at temperatures of kT ss J the effect of multi- 
partite entanglement cannot be neglected. In Section IV 
we demonstrate that our proofs also work for other mod- 
els by calculation similar thresholds for the XY model. 
Finally, we summarize our results and name some open 
problems. 
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II. NOTIONS OF MULTIPARTITE 
ENTANGLEMENT 

Let us start by clarifying the terms we will use to clas- 
sify multipartite entanglement in spin chains. For a pure 
state of a quantum systems of N parties we may ask 
for a given n < A, whether it is possible to cluster the 
parties into n groups, such that is a product state 
with respect to this partition, 

|V>> = I<M<H>I<£2>®-<8>K>. (1) 

If this is the case, we call the state l^) n-separable. A 
state which is A-separable is a product state with re- 
spect to all subsystems, these we call fully separable. If 
a state is not 2-separable, we call it genuine A -partite 
entangled. For a mixed state described by a density 
matrix g these terms can be extended through convex 
combination. If we can write g = ^2 i Pi\ipi){ipi\ with 
Pi > OiSi-Pi = 1 an d ^--separable \ipi) we call g n- 
separable. Physically, this means that a production of 
g requires only n-separable pure states and mixing [T^ . 

Besides asking for n-separability, we may also ask ques- 
tions like "Do two party entangled states suffice to create 
the mixed state g?" This leads to the following definition. 
We call a state producible by k -party entanglement (or 
fc-producible, for short) if we can write 

|V) = \<h) ® \<h) ® ... <S> \<j) m ) (2) 

where the \<pi) are states of maximally k parties. So, in 
this definition m > N/k has to hold. This definition 
means that it suffices to generate specific fc-party entan- 
glement to arrive at the state \ip). Conversely, we call a 
state containing genuine k -party entanglement if it is not 
producible by (fc — l)-party entanglement. This defini- 
tion can be extended to mixed states as before via convex 
combinations. Again, a mixed state which is ^-producible 
requires only the generation of /c-party entangled states 
and mixing for its production. Consequently, a mixed 
state g contains /c-party entanglement, iff the correla- 
tions cannot be explained by assuming the presence of 
(k — l)-party entanglement only. 

Obviously, there are some relations between the no- 
tions of fc-producibility and n-separability. For instance, 
the states containing A-party entanglement are just 
the genuine multipartite entangled states and the 1- 
producible states are the fully separable states. Further- 
more, a state of which some reduced state of m parties is 
genuine m-partite entangled, contains m-partite entan- 
glement, while the converse is in general not true |13| . 

In the thermodynamic limit Af> 1, however, the no- 
tion of n-separability for the total system becomes prob- 
lematic. This is mainly due to two reasons. First, it does 
not take into account how many particles are entangled. 
A biseparable state can be a product state for two parti- 
tions of equal size or just one qubit might be separated 
from a large genuine multipartite entangled state; these 
cases are not distinguished. Second, statements about 



n-separability require the exact knowledge of A, — a 
precondition which is usually not fulfilled in realistic sit- 
uations. In stark contrast, the notion of fc-producibility 
is designed to be sensitive to the question "How many 
qubits are entangled?" Also, deciding whether a state is 
fc-producible does not require the exact knowledge of N. 
If one controls only k <C A subsystems, one may still con- 
clude that the total state contains fc-party entanglement, 
e.g. if the reduced state is genuine fc-party entangled. 

III. THE HEISENBERG MODEL 

Let us start the discussion of spin models with the 
anti-ferromagnetic isotropic Heisenberg model. Here we 
assume a one-dimensional chain of A spin 1/2 particles 
with periodic boundary conditions. The Hamiltonian of 
this model is given by 

N 

Hh = j£oiM m) (3) 
i=i 

where a x denotes the Pauli matrix a x , acting on the i- 
th qubit and J > the coupling between the spins. For 
simplicity, we assume that A is even. For completeness, 
we always mention also known results for fully separable 
states [a, 13- We have the following theorems: 

Theorem 1. Let g be an A qubit state of a system 
described by the Heisenberg Hamiltonian in Eq. If 

(H H ) < -JN =: C m (4) 

there exists two neighboring qubits i,i + 1 in the chain 
such that the reduced state is a two-qubit entangled 
state. Furthermore, if 

(H H ) < -^-^JN « -1.618JA =: C R3 (5) 

then there exist three neighboring qubits i, i + l,i + 2 
such that the reduced state ft,i+i,i+2 of these qubits is 
genuine tripartite entangled. 

Proof. It suffices to prove these bounds for pure states. 
To prove first Eq. 10 we write Hh/J = W123 + W345 + 
1^567 + ... with W ijk = a^ ] a ( J ] + a { y l) 4 j) + af + 
o- ( x j) o- { x k) + 4 j) 4 fc) + ^Pa^ ■ The idea is to view the 
Wijk as entanglement witnesses, detecting genuine tri- 
partite entanglement on the qubits i,j,k. Let us denote 
Xi = (<Jx),XiXj = {(fx o~x ) etc. and consider a bisepa- 
rable state \ip) = \<j>i)\(f>jk) ■ We hav e \(Wjj k )\ = [ xj-Xj+yj- 

yj+Zi-Zj+XjXk+yjVk+ZjZkl < yxf+ y] + \xjXk+ 

yjVk + ZjZk\. For the two-qubit state \4>jk) the first term is 
the purity of the reduced state and invariant under local 
unitaries, while the second is maximal, when the 3x3 
correlation matrix A M „ = [ijV k for //, v = x, y, z is diago- 
nal. Expressing a general state in this form (see Eq. (14) 
in Ref. [14() leads to |(Wi,fc)| < 1 + v5. Since there are 
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N/2 of the Wijk, Eq. JSJ implies that one of the reduced 
three-qubit states cannot be biseparable. Finally, Eq. Q 
can be proved similarly, using \xt ■ Xj + j/j • yj + Zi ■ Zj | < 1 
for separable two-qubit states [a, la| ■ D 
Theorem 2. Let p be an N qubit state as in Theorem 
1. If g is 1-producible, then 

(H H ) > -JN =: C C 2 (6) 

holds, while for 2-producible states 

(H H ) > -\jN =: C C3 (7) 

holds. Thus, if (Hh) < Cc3 the state contains genuine 
tripartite entanglement. 

Proof. The bound Eq. ©.follows from Eq. @ and 
has already been derived in 0, . Let us first consider 
a two-producible pure state, where neighboring spins are 
allowed to be entangled. I.e., we have a state of the type 
\4') = |^i2>|<^34>"#Ar-i,iV>, where the state \<p 12 ) is a 
state of the qubits 1 and 2, etc. Let us define for the 
state \tp) two vectors with 6N real components each, via 

vr := ([1], [1:2], [2], [1], [5], [5:6], [6],..., [1]), 

v 2 := ([iV],[l],[3],[3:4],[4],[l],[7],...,[iV-l:7V]) 7 

where the symbols [...] stand always for three entries, 
namely [i] = Xi,yi,Zi, and = XiXj,yiyj, ZiZj and 

[1] = 1, 1, 1. Since the Vj have 6iV entries, in V\ (and 
similarly for v 2 ) some coefficients (like Hi) may appear 
twice, namely iff N € 4Z. It is straightforward to see that 
for the given state 2 • (Hh) = Jvi ■ v 2 holds. Now, we 
need X = rf+Vi+ z f + (xiX i+ i) 2 + (yiy l+ i) 2 + (ziz i+1 ) 2 + 
x i+i + Ui+i + z f+i ^ Xmax = 3, which holds for any two- 
qubit state g, due to (1 + x)/4 < Tr(g 2 ) < 1. Based on 
that, we have ||vi|| 2 < {N/2)(xmax + 3) = 3 AT. The same 
bound holds for ||i?2|| 2 - So, due to the Cauchy-Schwarz 
inequality we have \(H H )\ < l/2J||ui|| • ||t? 2 || < 3/2JN 
which proves the claim. 

Now we have to consider arbitrary two-producible 
states. A general two-producible pure state is always 
a tensor product of two-qubit states \4>ij) of the qubits i 
and j and single qubit states \4>k)- If for one of the two- 
qubit states \4>ij) the qubits i and j are not neighboring, 
we can replace it by two one-qubit reduced states gi (g> gj , 
since in this case, the Hamiltonian is only sensitive to 
the reduced states. Furthermore, if two neighboring sin- 
gle qubit states \<f>i) ® \4>j), appear, we replace them by 
\4>ij) 1 since \4>ij) is allowed to be separable. 

Thus it suffices to prove the bound for a state where 
entanglement is only present between neighboring qubits 
and the single qubit states are isolated in the sense that 
if \4>k) appears, then \4>k~i) and \4>k+i) do not appear, 
e.g. = \<pi2)\h)\M\<PM&)- Let M be the (even) 
number of isolated qubits. We can associate to any iso- 
lated qubit one neighboring two-qubit state (on the left 
or right) to form a three-party group. There is an am- 
biguity in doing that, and we can choose the three-qubit 



groups in such a way that the number of two-qubit groups 
between the three-qubit groups is always even. Then, 
the state can be viewed as a sequence gi, g 2 , 53, 34. •■ of 
M three-party groups and (AT— 3M)/2 two party groups, 
all in all there are (N — M)/2 groups. We can double this 
sequence by setting gr n-m) /2+k = 9k etc. In the example 
this grouping can be (12), (345), (678), (12), (345), (678). 
Now we define the vectors v\ and v 2 as follows: v\ col- 
lects terms from 31, 33, 35. ■■ while v 2 consists of terms 
from c/2, <?4, ■■■ In detail, we have 

vi = ([ffi], [51I52], [1], [52133], N, [53I54], [1], ...)> 
v 2 = ([1], [31I32], [32], [32I33], [1], [53 M, [34], ...), 

where [31] denotes the three terms XiXj,yiyj, ZiZj if 31 
is a two-qubit group and corresponds to the four terms 
XiXj,yiyj, ZiZj, (xjXk + yjyt + ZjZk) when 31 is a three- 
qubit group with the isolated qubit k. [31I32] etc. denotes 
the coupling terms between the groupings 31 and 32, in 
the Hamiltonian, e.g. Xi,yi,Zi for v\ and ^4.1,^+1,^+1 
for v 2 or vice versa. The symbol [1] denotes a sequence 
of three or four times "1", depending on whether in the 
other vector there is a two- or a three-qubit group. 

Again, we have 2 • (H) = Jv\ ■ v 2 and it remains to 
bound ||vi|| 2 . First, note that by construction V\ and v 2 
contain the same number of two and three-qubit groups, 
namely (N — 3M)/2 two-qubit groups and M three-qubit 
groups each. Also, both contain 3(A^ - 3M)/2 + AM 
times the "1". Then, note that for states of the type 
IVO = \<l>ij)\<j>k) theboundx 2 +y 2 + z 2 + (x l x j ) 2 + (y i y j ) 2 + 
(ziZj) 2 + (xj-Xk+yj-yk + Zj-z k ) 2 +xl+yl+zl < 5 is valid. 
So we have ||^|| 2 < (4 + 5)M + (3 + 3)(AT-3M)/2 = 3AT, 
which proves the claim. □ 

To start the discussion, first note that the bounds in 
Eqs. ijHl 0) are sharp. Equation J7J is saturated for the 
singlet chain | tf>) = \tp~)\tj}~)...\il>~) where = (|Q1) — 
\W))/V2. 

In general, any energy bound (Hh) > Cx corresponds 
to a certain temperature Tx • Below this temperature the 
state has with certainty some degree of entanglement. 
Corresponding to the Eqs. 0131 El 0) there are thus the 
temperatures Tr 2 ,Tr3,Tc 2 and Tc3 below which either 
reduced states of two or three parties are entangled or 
the total state contains two- or three-party entanglement. 
Obviously, T R2 = T C2 > T C 3 > T R3 has to hold here. 

Let us estimate these temperatures. For small N one 
can solve the Heisenberg model by diagonalizing Hh nu- 
merically |15| . Then, the threshold temperatures can di- 
rectly be computed. Results are shown in Fig. 1. As ex- 
pected, the values for Tci = Tr 2 coincide with the ones 
of Ref . |8j . The given values for Tc3 and Tr3 show that in 
the Heisenberg chain of ten spins at kT ~ J multipartite 
entanglement plays a role, namely at least one reduced 
state is genuine tripartite entangled and the total state 
contains tripartite entanglement. 

In the thermodynamic limit N 1 the ground state 
energy of the Heisenberg model is known to be E /N = 
-J(41n2 - 1) w -1.773J [Hj]. Thus three qubits can be 
found such that their reduced state is genuine tripartite 
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FIG. 1: Threshold temperatures Tr3,Tc2 and Tc3 for small 
Heisenberg chains up to ten qubits. See text for details. 



entangled. One can infer from numerical calculations [l5l 
IT?! ] that the threshold temperatures are determined by 
kTc3 « 1.61.7 and kT^ « 1.23J in agreement with the 
values of the ten-qubit Heisenberg chain. 

Finally, note that our results also shed light on the 
characterization of the ground state itself. For instance, 
they show that the ground state cannot be a GHZ state 
since for that state the reduced three-qubit states are sep- 
arable. Finally, note that when a state is translationally 
invariant, Eq. J^} guarantees that all reduced three-qubit 
states are genuine tripartite entangled. 



IV. THE XY MODEL 

Let us now investigate the isotropic XY model. The 
Hamiltonian of this model is given by 



HxY = jj2^ ) 4 +1) +-i; ) ^ +1) 



(8) 



Again, we assume periodic boundary conditions and an 
even number of spins. For this model we have: 

Theorem 3. Let g be an N qubit state whose dy- 
namics is governed by the Hamiltonian in Eq. ijSJl • If Q is 
one-producible, then 



(Hxy) > -JN 



(9) 



holds. If (Hxy) < — JN this implies that there are two 
neighboring qubits such that their reduced state is en- 
tangled. For two-producible states 



(Hxy) > -tjN 



(10) 



holds. If (Hxy) < —9/8JN the state contains thus tri- 
partite entanglement and if 

1 -\- f2 

(Hxy) < ^tt—JN « -1.207 JN (11) 



then there are three neighboring qubits such that their 
reduced state is genuine tripartite entangled. 

Proof. The proofs for the XY model are similar to the 
ones for the Heisenberg model, so we can make it short. 
Eq. © can be proved in the same manner as for the 
Heisenberg model [f|. Now, let us first show the bound 
(Unf for a state of the type \ip) — \<pi2)\<f>34,) —\<t>N-i,N) ■ 
Again, we define two vectors v\ and V2, with now AN 
entries via v x := ([1], [1:2], [2], [1], [5], [5:6], [6], [1], [1]) 
and v 2 := ([TV], [1], [3], [3:4], [4], [1], [7], [7:8], [N-hN]), 
where now [i] — Xi,yi and = af3i,2,ul3i,2 and [1] = 
1 1 with a — \/ 4/3 and [3ij = (xiXj + yiyj)/2. 



,a 

Again we have 2(Hxy) = Jv\ ■ V2 and due to the fact 



that xf + yf + 2a j! /3fj + yj + xj < £ max = 3 we have 
that ||vi|| 2 < N/2(£ max + 2/a 2 ) = 9/4N, which proves 
the claim. General two-producible states can be treated 
as shown in the proof of Theorem 2, now for states of the 
type \ipij)®\ipk) the bound xf +yf + 2a 2 /3f j +a 2 (x j -x k + 
yrVk) 2 +yl+xl < x 2 +yf+2a 2 (3l J +a 2 (x 2 j +y]) + l < 9/2 
must be used. Finally, Eq. © can be proved as Eq. 101, 

V2. 



using ,/x 2 + y 2 



< 1 



□ 



First note that again the bounds in Eqs. I|9I10|I are 
sharp. In Eq. I jlOjl equality holds for the state \ip) = 
\</>i2}\<t>3i)—-\(f>N-i,N) where \<j>i,i+i) (4>i,i+\ \ = 1/4- (I® 
1 - o- x ® u x - 1/2 • [o- y <g> a y + a z ® a z ] - \/3/2 ■ [1 ® 
Ox — o~x ® 1] ) for all odd i. This corresponds to a chain of 
non-maximally entangled two-qubit states. This also mo- 
tivates the definition of a and in the proof of Eq. i|10[l • 
a and were chosen such that for the state l^) we have 
vi = V2, thus the Cauchy-Schwarz inequality is sharp. 

In the thermodynamic limit the XY model is exactly 
solvable ^^]- The ground state energy is Eq/N = 
—AJ/tt w —1.273 J. Thus the ground state contains 
tripartite entanglement and fulfills the condition of 
Eq. 1)11(1 . Furthermore, the results of Ref. imply 
that if kT < 0.977J the mixed state contains tripartite 
entanglement and if kT < 0.668 J some reduced state is 
genuine tripartite entangled. 



V. CONCLUSION 

In conclusion, we showed for two important spin mod- 
els that the internal energy can be used as a signature for 
the presence of multipartite entanglement in these mod- 
els. Based on this, we computed threshold temperatures 
below which any realistic description of the system can- 
not neglect the effects of multipartite entanglement. Our 
results may stimulate the research on entanglement in 
phase transitions, since they suggests the use of multipar- 
tite entanglement measures as a tool for the investigation 
of phase transitions in these regimes |19|. 

A natural continuation of the present work is the ex- 
tension of the Theorems presented here to other spin 
systems, e.g. higher dimensional or frustrated systems. 
Furthermore, it is very tempting to analyze the results of 
performed experiments whether they can be interpreted 
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as giving evidence for multipartite entanglement, simi- 
larly as it has been done for bipartite entanglement in 
Ref. [2D|. Here, we leave this as an open problem. 
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